Doping dependence of Meissner effect in cuprate superconductors 
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Within the t-t'-J model, the doping dependence of the Meissner effect in cuprate superconductors 
is studied based on the kinetic energy driven superconducting mechanism. Following the linear 
response theory, it is shown that the electromagnetic response consists of two parts, the diamagnetic 
current and the paramagnetic current, which exactly cancels the diamagnetic term in the normal 
state, and then the Meissner effect is obtained for all the temperature T < T c throughout the 
superconducting dome. By considering the two-dimensional geometry of cuprate superconductors 
within the specular reflection model, the main features of the doping and temperature dependence 
of the local magnetic field profile, the magnetic field penetration depth, and the superfluid density 
observed on cuprate superconductors are well reproduced. In particular, it is shown that in analogy 
to the domelike shape of the doping dependent superconducting transition temperature, the maximal 
superfluid density occurs around the critical doping 8 ~ 0.195, and then decreases in both lower 
doped and higher doped regimes. 

PACS numbers: 74.25.Ha, 74.25.Nf, 74.20.Mn 



I. INTRODUCTION 

One of the characteristic features of superconductors 
is the so-called Meissner effect^, i.e., a superconductor 
is placed in an external magnetic field B smaller than 
the upper critical field B c , the magnetic field B pene- 
trates only to a penetration depth A (few hundred nm 
for cuprate superconductors at zero-temperature) and is 
excluded from the main body of the system. This mag- 
netic field penetration depth is a fundamental parameter 
of superconductors, and is closely related to the super- 
fluid density p^&, which is proportional to the squared 
amplitude of the macroscopic wave function, and there- 
fore describes the superconducting (SC) charge carriers. 
Furthermore, the variation of the magnetic field pene- 
tration depth (then the superfluid density) as a function 
of doping and temperature provides information crucial 
to understanding the details of the SC stated. In par- 
ticular, since the compounds of cuprate superconductors 
are doped Mott insulators with the strong short-range 
antiferromagnetic correlation dominating the entire SC 
phase 4 , the magnetic field can be also used to probe the 
doping and momentum dependence of the SC gap and 
spin structure of the Cooper pair—. This is why the first 
evidence of the d-wave Cooper pairing state in cuprate 
superconductors was obtained from the earlier experi- 
mental measurement for the magnetic field penetration 
depths. Since then this d-wave SC state remains one of 
the cornerstones of our understanding of the physics in 
cuprate superconductors^. 

Experimentally, by virtue of systematic studies using 
the muon-spin-rotation measurement technique, some es- 
sential features of the evolution of the magnetic field pen- 
etration depth and superfluid density in cuprate super- 
conductors with doping and temperature have been es- 
tablished now for all the temperature T < T c throughout 
the SC dome: (1) the magnetic field screening is found 
to be of exponential character^, in support of a local 
(London- type) nature of the electrodynamic response 2 ; 



(2) the magnetic field penetration depth is a linear tem- 
perature dependence at low temperatures except for 
the extremely low temperatures where a strong devia- 
tion from the linear characteristics (a nonlinear effect) 
appears^ - — ; (3) the doping dependence of the super- 
fluid density is observecU 2 - - — , where the superfluid den- 
sity exhibits a peak around the critical doping 8 ~ 0.19, 
and then decreases at both lower doped and higher doped 
regimes. This in turn gives rise to the domelike shape 
of the doping dependence of the SC transition tempera- 
ture. In particular, it has been shown experimentally 13 
that the peak of the superfluid density at the critical 
doping 8 ~ 0.19 is a common feature of the hole-doped 
cuprate superconductors. Theoretically, to the best of 
our knowledge, all theoretical calculations of the exper- 
imental measurements of the magnetic field penetration 
depth and the related superfluid density in cuprate su- 
perconductors performed so far are based on the phe- 
nomenological d-wave Bardeen-Cooper-Schrieffer (BCS) 
formalism 1 ^—. In the local limit, where the magnetic 
field penetration depth A is much larger than the coher- 
ence length C, i.e., A 3> it has been shown&i£ that 
the simple d-wave pairing state gives the linear tempera- 
ture dependence of the magnetic field penetration depth 
AA(T) = A(T) - A(0) oc T/A at low temperatures, with 
Aq is the zero-temperature value of the d-wave gap am- 
plitude. However, the characteristic feature of the d-wave 
energy gap is the existence of the gap nodes on the Fermi 
surface, which can lead to the nonlinear effect of field on 
the penetration depth (then superfluid density) at the 
extremely low temperatures, associated with a field in- 
duced increase in the density of available excitation states 
located around the gap nodes on the Fermi surface^— . 
This follows from a fact that the nonlocal effect is closely 
related to the divergence of the coherence length ( at 
the gap nodes on the Fermi surface, as the coherence 
length £ varies in inverse proportion to the value of the 
energy gap, then this nonlocal effect at the extremely low 
temperatures can lead to a nonlinear temperature depen- 
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dence of the magnetic field penetration depth in the clean 
limiti^o. 

In our recent workSi based on the nearest neighbors 
hopping t-J model, the weak electromagnetic response 
in cuprate superconductors has been discussed within the 
kinetic energy driven SC mechanism 22 . However, it has 
been shown experimentally^ 2 ^ that although the highest 
energy filled electron band is well described by the near- 
est neighbors hopping t-J model, the overall dispersion 
may be properly accounted by generalizing the nearest 
neighbors hopping t-J model to include the second- and 
third- nearest neighbors hopping terms t' and t". Further- 
more, the experimental analysis 2 ^ showed that the SC 
transition temperature (then the superfluid density) for 
different families of cuprate superconductors is strongly 
correlated with t' . In this paper based on the t-t'-J 
model, we study the doping dependence of the Meissner 
effect in cuprate superconductors for all the temperature 
T < T c throughout the SC dome, where one of our main 
results is that the superfluid density increases with in- 
creasing doping in the lower doped regime, and reaches 
a maximum around the critical doping 8 ~ 0.195, then 
decreases in the higher doped regime. 

The rest of this paper is organized as follows. The ba- 
sic formalism is presented in Section [XXI where within the 
t-t'-J model, we evaluate the diamagnetic and paramag- 
netic components of the response kernel function based 
on the kinetic energy driven SC mechanism, and then 
obtain explicitly the Meissner effect in cuprate supercon- 
ductors for all the temperature T < T c throughout the 
SC dome. Within this theoretical framework, we discuss 
the basic behavior of cuprate superconductors in a weak 
electromagnetic field in Section IIII1 and qualitatively re- 
produce all the main features of the doping and tem- 
perature dependence of the local magnetic field profile, 
the magnetic field penetration depth, and the superfluid 
density. Finally, we give a summary and discussions in 
Section [IV] In Appendix [X] we give the explicit forms of 
the paramagnetic part of the response kernel function at 
the temperatures T — and T = T c , respectively. 



II. THEORETICAL FRAMEWORK 

In cuprate superconductors, the characteristic feature 
is the presence of the Cu02 planed. It has been argued 
that the essential physics of the doped Cu02 plane^^ 
is properly accounted by the t-t'-J model on a square 
lattice. However, for discussions of the doping and tem- 
perature dependence of the Meissner effect in cuprate 
superconductors, the t-t'-J model can be extended by 
including the exponential Peierls factors as, 

Irja la 
If]' a Ifj 



supplemented by an important on-site local constraint 
J2 a C\ a Ci a < 1 to remove the double occupancy, where 
fj = ±x, ±y, fj' = ±x ± y, C\ a (Cia) is the electron cre- 
ation (annihilation) operator, Si = (Sf , Sf , Sf) are spin 
operators, and /i is the chemical potential. The expo- 
nential Peierls factors account for the coupling of elec- 
trons to the weak external magnetic fiel d 26 ' 27 in terms 
of the vector potential A(Z). To incorporate the electron 
single occupancy local constraint, the charge-spin sepa- 
ration (CSS) fermion-spin theory 2 ^ has been proposed, 
where the physics of no double occupancy is taken into 
account by representing the electron as a composite ob- 
ject created by Cif — h\^Sf and = hj^Sf , with the 
spinful fermion operator hi a — e~ l ^ la hi that describes 
the charge degree of freedom of the electron together 
with some effects of spin configuration rearrangements 
due to the presence of the doped hole itself (charge car- 
rier), while the spin operator Si represents the spin de- 
gree of freedom of the electron, then the electron single 
occupancy local constraint is satisfied in analytical calcu- 
lations. In particular, it has been shown that under the 
decoupling scheme, this CSS fermion-spin representation 
is a natural representation of the constrained electron 
defined in the Hilbert subspace without double electron 
occupancy 2 ^. In this CSS fermion-spin representation, 
the t-t'-J model ([T]) can be expressed as, 

H = t^e-^^HhU^S+S^ 
if) 

+ h l+fj'i, h ii s r s kv>) 

- n h'fohfo + Jeff S; • Si +fj , (2) 

la If) 

where J c g = (1 — 8) 2 J, and 8 — (h' llT hi a ) = (h]hi) is the 
charge carrier doping concentration. 

As in the conventional superconductors, the key phe- 
nomenon occurring in cuprate superconductors in the SC 
state is the pairing of charge carriers^. The system 
of charge carriers forms pairs of bound charge carriers 
in the SC state, while the pairing means that there is 
an attraction between charge carriers. What is the ori- 
gin of such an attractive force? Recently, the kinetic 
energy driven SC mechanism has been developed based 
on the t-J model 2 ^, where the charge carrier-spin inter- 
action from the kinetic energy term in the t-J model 
((2|) induces a d-wave charge carrier pairing state by ex- 
changing spin excitations in the higher power of the dop- 
ing concentration, then the SC transition temperature 
is identical to the charge carrier pair transition temper- 
ature. Furthermore, this SC state is the conventional 
BCS-like with the d-wave symmetr y 29 ! 30 , so that the ba- 
sic d-wave BCS formalism is still valid in quantitatively 
reproducing all main low energy features of the SC co- 
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herence of the quasiparticle peaks in cuprate supercon- 
ductors, although the pairing mechanism is driven by the 
kinetic energy by exchanging spin excitations. Following 
these previous discussion s 22 ! 29 ' 30 , the full charge carrier 
Green's function in the zero magnetic field case can be 
obtained explicitly in the Nambu representation as, 



G(k,iw„) = Z hF 



w„t + £ k r 3 - A h z(k)n 



(3) 



where To is the unit matrix, t\ and T3 are Pauli matri- 
ces, the renormalized charge carrier excitation spectrum 
£k = ZhFCk, with the mean- field charge carrier excita- 
tion spectrum £k = Ztxijk — Zt'xil'^ ~ Mi the spin 
correlation functions xi = {S^S~ + ~), \2 = (S^S^ +fjl ), 

7k = (1/^)E^^, 7k = (WEfl/e*^, Z is the 
number of the nearest neighbor or second-nearest neigh- 
bor sites, the renormalized charge carrier d-wave pair 
gap Ahz(k) = Zhp Ah (k), where the effective charge 
carrier d-wave pair gap Ah(k) = Ah(cosfc x — cosfc y )/2, 
and the charge carrier quasiparticle spectrum Ehk = 



|Ahz(k)| 2 , while the effective charge carrier gap 



parameter Ah and the quasiparticle coherent weight ZhF 
have been determined self-consistently along with other 
equation s 29 ' 30 , and then all order parameters and chemi- 
cal potential have been determined by the self-consistent 
calculation. 



A. Linear response approach in the presence of a 
weak external magnetic field 

In cuprate superconductors, an external magnetic field 
generally represents a large perturbation on the system, 
then the induced field arising from the superconductor 
cancels this external magnetic field over most of the sys- 
tem. In this case, the net field acts only near the sur- 
face on a scale of the magnetic field penetration depth, 
and then it can be treated as a weak perturbation on 
the system as a whole 2 . This is why the Meissner ef- 
fect can be successfully studied within the linear response 
approac h 31 ! 32 , where the linear response current density 
J M and the vector potential A v are related by a nonlocal 
kernel of the response function K^u as, 



(4) 



with the Greek indices label the axes of the Cartesian 
coordinate system. The kernel of the response function 
in Eq. (Q| plays a crucial role for the discussion of the 
doping dependence of the Meissner effect in cuprate su- 
perconductors, and can be separated into two parts as, 



(5) 



where K^J and Kjfi) are the corresponding diamagnetic 
and paramagnetic parts, respectively, and are closely re- 



lated to the electron current density in the presence of 
the vector potential A u . 

The vector potential A (then the weak external mag- 
netic field B — rot A) has been coupled to the elec- 
trons, which are now represented by Cj-p = hj^S^ and 

Cix = hnSf hi the CSS fermion-spin representation. In 
this case, the electron current operator can be obtained 
in terms of the electron polarization operator, which is 
a summation over all the particles and their positions, 
and can be expressed explicitly in the CSS fermion-spin 
representation as, 



P = -e J2 RiCfa-CW = e J2 R A^ 



(6) 



then the electron current operator is obtained by evalu- 
ating the time-derivative of this polarization operator ([6]) 



iet 



in 



+ h n h\ +fii S l S+ + .) 



iet' 



If)' 



+ h n h\ +f) , i S l 5+ ry ) 



(7) 



In the linear response approach, this electron current op- 
erator is reduced as j = , with the corresponding 
diamagnetic (d) and paramagnetic (p) components of the 
electron current operator are given by, 



j (d) = iE«A(i)'«KW s i 



r, 2 



in 



+ h n h l +fl i s i 5 i+r» 



eH' 
h 2 



+ h n h l+f,'i s i s Ufy) > 



(8a) 



;(p) 



iet 



in 



If)' 



(8b) 



respectively. Obviously, the diamagnetic component of 
the electron current operator in Eq. ([8a]) is proportional 
to the vector potential, and in this case, we can obtain 
the diamagnetic part of the response kernel directly as, 

K$(<i,u) = -^-(xi<M-2 X2 2 i')<V = -^-6^, (9) 
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where the charge carrier particle-hole parameters <f>\ = 
(h\ a h i+fia ) and (f> 2 = (hl a h l+ ,y a }, while A^ 2 = 
—4e 2 (xi4>it — 2x2<t>2t / )/h 2 is the London penetration 
depth, and is doping and temperature dependent. 

However, the paramagnetic part of the response kernel 
is more complicated to calculate, since it involves evalua- 
tion of the following electron current-current correlation 
function (polarization bubble), 

p»M,t) = -(r r {jW(q.T-)^ } (-q»o)}) J (10) 

then the paramagnetic part of the response kernel 
K/fi) (q, uj) can be obtained as K$ (q, uj) = -P A ,„(q, w). In 
the CSS fcrmion-spin approach, the paramagnetic com- 
ponent of the electron current operator in Eq. (|8bl) can 



be decoupled as, 

J - ^ / j 'I'H+fKj'Ha i- ^ 2-^ ' l+fi'<r la 

If/a Ifj' a 

+ l -^f^{ s ts W r + s^st +il ), (ii) 

where the third and fourth terms in the right-hand side 
refer to the contribution from the electron spin, and can 
be expressed explicitly as, 



f ^ ^[(^t^i+p, + &i Sj+p) (Si'Si-p, + &i Slip)] — ' ft[(S^_ a S l + S l+fi Si~) - (s^s l+jl + S t S^_p)} 



l/'l 



o, 



Y^[(x + y)(S+S l+£+ ~ + S t Sf + - + ~) -(x + y^SfS^^. + S t S+_ & _-) + (4 - y)(S+ S l+i _$ + S t Sf +£ _.) 



In this case, the majority contribution for the paramag- 
netic component of the electron current operator comes 
from the electron charge, and then the paramagnetic 
component of the electron current operator in Eq. (|11[) 
can be expressed explicitly as, 

la 

- hl^hi+fi+ta + hi+ji-tiv)]. (12) 

For the convenience in the following discussions, the para- 
magnetic component of the electron current operator can 



be rewritten in the Nambu representation in terms of 
the charge carrier Nambu operators ^ = (h^, h-kij 

and ^k+q = ^k+qt: ^-k-qi) • Moreover, since the 
density operator is summed over the position of all 
particles, then its Fourier transform can be obtained 

as p(q) = (e/2)£ k< X<A+ qCT = ( e / 2 )E k ^3*k +q . 
In this Nambu representation, the paramagnetic four- 
current operator can be represented as, 

^(qH^E^k.k + q^k+q. ( 13 ) 
where the bare current vertex, 



-|e2*{sin(fe M + 2x2*' E cos (| ft,) cos(fc„ + \q v )] 

7 M (k + q,k) = { -i(2xit') cos(fc M + ±ft,) £ sin ft, s\n(k v + \q v )}T Q for /i ^ 0, 

|r 3 for /i = 0. 



(14) 



As in the previous discussions—, we are calculating 
the polarization bubble with the paramagnetic current 
operator (p~3|) . i.e., bare current vertices (fM)) . but Green 
function ([3]). As a consequence, we do not take into ac- 
count longitudinal excitations properly^!, the obtained 



results are valid only in the gauge, where the vector po- 
tential is purely transverse, e.g. in the Coulomb gauge. 
In this case, the correlation function ([TO]) can be obtained 
in the Nambu representation as, 



P M „(q,iw„) = ^^7 M (k + q,k)7*(k + q,k)i^Tr[G(k + q,iw n + zz/ m )G(k,w m )] . (15) 

Substituting the Green's function ([3]) into Eq. (|15l) . we then obtain the paramagnetic part of the response kernel in 
the static limit (w ~ 0) as, 

K${<\^) = ^^7/,(k + q,k)7:(k + q,k)[£ 1 (k,q)+i 2 (k,q)] = 4P)(q,0)V, (16) 

k 

with the functions Za(k, q, ui) and ^2(k, q, ui) are given by, 

j ri x 7 2 (, , gk+qlk + A hz (k + q)A h z(k) \ n F (E hk ) - n F (£ hk+q ) 

Ll\H-,<l) — ^hF 1 H 5 — ^ 5 5 ' l i ' a v 

\ -C/hk-C'hk+q / -C/hk — -C/hk+q 

L 2 (k,q) = Z 2 F ( 1 - ^Wk + A hz (k + q)A hz (k)\ n F (E hk )+n F (E hk+cl )-^ 

\ EhkEhk+q ) £hk + -Ehk+q 



respectively. In this case, the kernel of the response func 
tion in Eq. ([5]) is now obtained from Eqs. ((9]) and ((TBI 

as, 



A^(q,0) = 



1 



4 P M } (q,o) 



(18) 



It should be emphasized that in the present CSS fermion- 
spin theor y 28 ' 29 , these charge carrier h\ and spin S/ are 
gauge invariant, and in this sense they are real and can 
be interpreted as physical excitations^. Furthermore, as 
shown in Eq. © and Eq. fT2")l , the electron polarization 
operator and the related electron current operator are 
identified with the corresponding charge carrier polariza- 
tion operator and charge carrier current operator, since 



the electron single occupancy local constraint is satisfied 
in the CSS fermion-spin approach. 



Doping dependence of the Meissner effect in 
the long wavelength limit 



B. 



With the help of the response kernel function (|18l) . 
we now discuss the doping and temperature dependence 
of the Meissner effect in cuprate superconductors. In 
particular, in the long wavelength limit, i.e., |q| — > 0, the 
function ^(k, q — > 0) vanishes, then the paramagnetic 
part of the response kernel can be obtained explicitly as, 



A(P)(q^0,0) 



4 e 2 i 

2Z hF-^2~^X! Sin2 MXl* - 2 X2t' COS k x 



I 2 lim 

q^O 



"F(-Ehk) " n F (E hkA 
£?hk — -Ehk+q 



(19) 



At zero-temperature T = 0, Kyy (q —> 0, 0)|t=o = (see 
Appendix |A]) . In this case, the long wavelength electro- 
magnetic response is determined by the diamagnetic part 
of the kernel only. On the other hand, at the SC transi- 



tion temperature T = T c , the effective charge carrier gap 
parameter Ah|r=T c = 0. In this case, the paramagnetic 
part of the response kernel in the long wavelength limit 
can be evaluated as (see Appendix |A"| . 



A(P)(q^0,0)| T=Tc = 



2Z hF-TT T? 2^ sm k v \* lt ~ 2 * 2t cos k x Y lim = = ^ 

" iV k Sk - ?k+q 



"Al' 



(20) 
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which exactly cancels the diamagnetic part of the re- 
sponse kernel (|9|) , and then the Meissner effect in cuprate 
superconductors is obtained for all T <T C . To show this 
point clearly, we define the effective superfluid density 
n s (T) at temperature T in terms of the paramagnetic 
part of the response kernel as, 



ifCp)(q^0,0) = -i 



1 - 



n«(T) 
n s (0) 



(21) 



and then the kernel of the response function in Eq. (TT5)) 
can be rewritten as, 



i^„(q-K),0) = 



f 



n s (T) s 
n s (0) V ' 



(22) 



where the ratio n s (T)/n s (0) of the effective superfluid 
densities at temperature T and zero-temperature is given 
by, 



n s {T) 
n a (0) 



f - 2A|Z 2 F ^^ J2 sin2 k vlxit - 2 X 2t' cos k x f 



fie 



( e /3E hk + 1)2 



In cuprate superconductors, although the values of J, t, 
and t' are believed to vary somewhat from compound 
to compound^, however, as a qualitative discussion, the 
commonly used parameters in this paper are chosen as 
t/J = 2.5, t'/t = 0.3, and J = I000K. In this case, 
we plot the effective superfluid density n s (T)/n s (0) as a 
function of temperature T for the doping concentration 
6 = 0.09 (solid line), 6 = 0.12 (dashed line), and S = 0.15 
(dash-dotted line) in Fig. [TJ Obviously, it is shown that 
the effective superfluid density diminishes with increas- 
ing temperatures, and disappears at the SC transition 
temperature T c , then all the charge carriers are in the 
normal fluid for the temperature T > T c . In summary, 
we have found the following results within the kinetic en- 
ergy driven SC mechanism: (1) the doping dependence 
of the Meissner effect in cuprate superconductors is ob- 
tained for all temperature T < T c throughout the SC 
dome; (2) the electromagnetic response kernel goes to 
the London form in the long wavelength limit [see Eq. 
([22])] ; (3) although the electromagnetic response kernel 
is not manifestly gauge invariant within the present bare 
current vertex (TH)) . it has been shown that the gauge 
invariance is kept within the dressed current vertes^i. 



III. QUANTITATIVE CHARACTERISTICS OF 
THE MEISSNER EFFECT 



In this section, we discuss the basic behavior of cuprate 
superconductors in the presence of a weak electromag- 
netic field. As seen from Eq. (U), once the response ker- 
nel is known, the effect of a weak electromagnetic 
field can be quantitatively characterized by experimen- 
tally measurable quantities such as the local magnetic 
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FIG. 1: The effective superfluid density as a function of tem- 
perature T for the doping concentration 8 = 0.09 (solid line), 
S = 0.12 (dashed line), and 5 = 0.15 (dash-dotted line) with 
parameters t/J = 2.5, t'/t = 0.3, and J = 1000K. 
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field profile, the magnetic field penetration depth, and 
the superfluid density. However, the result we have ob- 
tained the response kernel in Eq. (fT5|) [then the effec- 
tive superfluid density (|2"3")) ] can not be used for a direct 
comparison with the corresponding experimental data of 
cuprate superconductors because the kernel function de- 
rived within the linear response theory describes the re- 
sponse of an infinite system, whereas in the problem of 
the penetration of the field and the system has a surface, 
i.e., it occupies a half-space x > 0. In such problems, 
it is necessary to impose boundary conditions for charge 
carriers. This can be done within the simplest specular 
reflection mode l 35 ' 36 with a two-dimensional geometry of 
the SC plane. In this case, one may use the response 
kernel for an infinite space, however, it is necessary 
to extend the vector potential A(r) in an even manner 
through the boundary. If the externally applied magnetic 
field is perpendicular to the ab plane, we may choose 
Ay(x) along the y axis. Following the Maxwell equation, 

rotB = rot rot A = graddivA - V 2 A = ^ J, (24) 

it is shown clearly that the extension of the vector poten- 
tial in an even manner through the boundary implies a 
kink in the A y (x) curve. In other words, if the externally 
applied magnetic field B is given at the system surface, 
i.e., (dA y (x)/dx)\ x=+ Q = B, while (dA y (x)/dx)\ x= - = 
—B, this leads to a fact^ 5 - that the second derivative 
(d 2 A y (x)/d 2 x) acquires a correction 2B5(x), i.e., 



d 2 A y (x) 
d 2 x 



2BS(x) — hqJ v 



(25) 



where the transverse gauge div A = has been adopted. 
In the momentum space, this equation (1251) can be ex- 
pressed as q 2 A y (q) = fj,oJ y (q) — 2B. Substituting this 
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Fourier transform of Eq. (1251) into Eq. ((4]), and solving 
for the vector potential we obtain, 



Mi) 



-IB 



fi K yy (ci) +ql 



(26) 



Since the vector potential has only the y component, the 
non-zero component of the local magnetic field h = rot A 
is that along the z axis and h z (q) = iq x A y (q). With 
the help of Eq. ([2T5)) . we can obtain explicitly the local 
magnetic field profile as, 



h z (x) 



D 



q x sin(g x a;) 



fJ, K yy (q x , 0,0) +ql' 



(27) 



which therefore reflects the measurably electromagnetic 
response in cuprate superconductors. For the conve- 
nience in the following discussions, we introduce a charac- 
teristic length scale ao = \J ti 1 a / noe 2 J . Using the lattice 
parameter a ~ 0.383nm for the cuprate superconduc- 
tor YBa2Cu307_ a , this characteristic length is obtain as 
a w 97.8nm. In Fig. [21 we plot the local magnetic 
field profile (|2"?| as a function of the distance from the 
surface at temperature T = 2K for the doping concen- 
tration S — 0.09 (solid line), 5 = 0.12 (dashed line), and 
S = 0.15 (dash-dotted line) in comparison with the corre- 
sponding experimental result^ of the local magnetic field 
profiles for the high quality YBa2Cu307_ y (inset). If an 
external field B — 8.82mT is applied to the system just 
as it has been done in the experimental measurement 8 , 
then the experimental result for YBa2Cu307_ a is well 
reproduced. In particular, our theoretical results per- 
fectly follow an exponential field decay as expected for 
the local electrodynamic response. This is different from 
the conventional superconductors, where the local mag- 
netic field profile in the Meissner state shows a clear de- 
viation from the exponential field decay^. The exponen- 
tial character of the local magnetic field profile has been 
observed experimentally on different families of cuprate 
superconductors 7 -^ , in support of a local (London- type) 
nature of the electrodynamics 2 . 



A. Doping and temperature dependence of the 
magnetic field penetration depth 

In the following discussions, we do not analyze the be- 
havior of the filed in the general case, and only discuss 
the doping and temperature dependence of the magnetic 
field in-plane penetration depth A(T) and the related in- 
plane superfluid density p s (T). The magnetic field in- 
plane penetration depth is defined in terms of the local 
magnetic field profile (|27)) as, 



A(T) = i = t? 

B J 7T J ^ o K yy (q x ,0, 



yy(qx,0,0) + q2' 



(28) 



At zero-temperature, the calculated magnetic field pen- 
etration depths are A(0) « 239.17nm, A(0) w 234.76nm, 
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FIG. 2: The local magnetic field profile as a function of the 
distance from the surface at temperature T = 2K for the 
doping concentration S = 0.09 (solid line), S = 0.12 (dashed 
line), and 5 — 0.15 (dash-dotted line) with parameters t/J = 
2.5, t'/t = 0.3, and J = 1000K. Inset: the corresponding 
experimental result for YBa20u3O7_ H taken from Ref. [g. 



and A(0) « 224.44nm for the doping concentrations 
S = 0.14, 5 = 0.15, and S = 0.18, respectively, which 
are consistent with the values of the magnetic field pen- 
etration depth A s» 156nm ~ 400nm observed for differ- 
ent families of cuprate superconductors in different dop- 
ing concentrations^ ^ 14 ' 37 . On the other hand, at the 
SC transition temperature T = T c , the kernel of the re- 
sponse function lf M „(q — > 0, 0)|t=t c = 0. In this case, 
we obtain the magnetic field penetration depth from Eq. 
(|2"5)) as A(T C ) = oo, which reflects that in the normal 
state, the external magnetic field can penetrate through 
the main body of the system, therefore there is no the 
Meissner effect in the normal state. For a better un- 
derstanding of the unusual behavior of the temperature 
dependence of the magnetic field penetration depth A(T), 
AA(T) = A(T) — A(0) as a function of temperature T for 
the doping concentration 5 = 0.14 (solid line), 5 = 0.15 
(dashed line), and S — 0.18 (dash-dotted line) is plot- 
ted in Fig. [3] in comparison with the corresponding ex- 
perimental results^ of YB^CuaC^-y (inset). Our result 
shows clearly that in low temperature, the magnetic field 
penetration depth AA(T) exhibits a linear temperature 
dependence, however, it crosses over to a nonlinear be- 
havior in the extremely low temperatures, in good agree- 
ment with experimental observation in nominally clean 
crystals of cuprate superconductors^ 7 - - — . In compari- 
son with our previous discussions 21 , our present results 
also show that the good agreement can be reached by in- 
troducing the second-nearest neighbors hopping t' in the 
nearest neighbors hopping t-J model. It should be em- 
phasized that the present result for cuprate superconduc- 
tors is much different from that in the conventional super- 
conductors, where the characteristic feature is the exis- 
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tence of the isotropic energy gap A s , and then AA(T) ex- 
hibits an exponential behavior as AA(T) oc exp(— A s /T). 
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FIG. 3: Temperature dependence of the magnetic field pen- 
etration depth AA(T) for the doping concentration S — 0.14 
(solid line), S — 0.15 (dashed line), and 8 — 0.18 (dash- 
dotted line) with parameters t/J — 2.5, t'/t — 0.3, and 
J — 1000K. Inset: the corresponding experimental data for 
YBa2Cu307_i, taken from Ref. [{|. 



B. Doping and temperature dependence of the 
in-plane superfluid density 

Now we turn to discuss the doping and temperature de- 
pendence of the in-plane superfluid density p s (T) , which 
is a measure of the phase stiffness, and can be obtained 
in terms of the magnetic field in-plane penetration depth 
A(T) as, 



Ps(T) 



1 



A 2 (T)- 



(29) 



In this case, we have performed firstly a calculation for 
the doping dependence of the zero-temperature super- 
fluid density p s (0) for all levels of doping, and the result 
is plotted in Fig. 0]in comparison with the corresponding 
experimental data^ for Yo.8Cao.2Ba2(Cui_ z Zn 2 )307_,5 



and Tli-yPbyS^Cai-zY-rC^Or (inset). It is shown 
clearly that the zero-temperature superfluid density in- 
creases with increasing doping in the lower doped regime, 
and reaches a maximum (a peak) around the critical 
doping S « 0.195, then decreases in the higher doped 
regime, in good agreement with the experimental results 
of cuprate superconductors^^— . Since the superfluid 
density p s (T) [then the magnetic field penetration depth 
A(T)] is related to the current-current correlation func- 
tion, the charge carrier pair gap parameter is relevant as 
shown in Eqs. (fTrJ)) . (|T7|) . and (fP5)) , i.e., the variation 
of the superfluid density with doping and temperature 
is coupled to the doping and temperature dependence of 
the charge carrier pair gap parameter A^z in cuprate su- 
perconductors. In this case, our present domelike shape 
of the doping dependent superfluid density also is a nat- 
ural consequence of the domelike shape of the doping 
dependence of the SC transition temperature (then the 
charge carrier pair gap parameter) in the framework of 
the kinetic energy driven SC mechanism^ 2 ., where the 
the maximal SC transition temperature (then the charge 
carrier pair gap parameter) occurs in the optimal dop- 
ing, and then decreases in both underdoped and over- 
doped regimes. However, the calculated SC transition 
temperature T c [then the zero-temperature gap param- 
eter Ah(0)] exhibits the maximal value at the optimal 
doping ^optimal ~ 0.1 52S, therefore there is a difference 
between the optimal doping ^optimal ~ 0.15 [the maxi- 
mal Ah(0) value] and the critical doping Critical ~ 0.195 
[the highest p s (0) value]. This difference can be un- 
derstood within the present theoretical framework. In 
the domelike shape of the doping dependence of the gap 
parameter Ah(0), the gap parameter Ah(0) reaches its 
maximal value at the optimal doping <5 ptimai ~ 0.15, 
where the doping-derivative of Ah(0) can be obtained as 
(dA h (0)/d<5)|«5 = 5 optimal = 0. On the other hand, at the 
critical doping (5 cr iti ca i ~ 0.195, the peak of the super- 
fluid density /0 S (0) appears, where the doping-derivative 
of p s (0) can be obtained as {dp s (0)/dS)\s=s CIitical = 
0. According to the definition of the superfluid den- 
sity p s (T) in Eq. Jgg, (dp s (0)/d<5)| 5=5critical = is 
equivalent to (dA(0)/d5)| < 5 = 5 crWca] = 0. In this case, 
(dA(0)/d<5)|«5 = 5 critical = can be expressed from Eq. ([28]) 
as, 



"dA(0)" 



2p Q 



dq x 



dK yy (q x ,0,0)~ 



s=s c . 



[MoK yy (q x ,0,0) + q 2 



(30) 



<5=<5critical 



then with the help of Eqs. ©, ([16). and (JT7J), it is 
straightforward to find that when {dp s (0)/d8)\s=s„ itical = 
0, (dAh(0)/d<5)|a = 5 critical ^ 0, since the diamagnetic part 



of the response kernel in Eq. the spin correla- 

tion functions xi = ( S i~ S i+f,) and X2 = {SfS^ + ~,), 
and the charge carrier particle-hole parameters <fi\ — 



9 



(h'^hi+fjcr) and <f) 2 = (hl a hi+^ a ) are also doping de- 
pendent. This leads to the difference between the op- 
timal doping <5 ptimai ~ 0.15 for the zero-temperature 
gap parameter (then the SC transition temperature) 
and the critical doping <5 C riticai ~ 0.195 for the zero- 
temperature superfluid density. In particular, it is found 
that (dAh(0)/d5)|,5 =( 5 critical < 0, indicating that the criti- 
cal doping locates at the slightly overdoped regime. Fur- 
thermore, it should be emphasized that the early ex- 
perimental data observed from cuprate superconductors 
show that the superfluid density p s (0) in the underdoped 
regime vanishes more or less linearly with decreasing 
doping concentration^. Later, a clear deviation from 
this linear relation between the superfluid density p s {0) 
and doping concentration has been observed in the un- 
derdoped regimeii~— . In particular, the recent exper- 
imental measurement^ on the cuprate superconductor 
YBa2Cu307_ a indicate that the superfluid density p s (0) 
is, in actual fact, linearly proportional to the doping con- 
centration in the low doping range (S « 0.054 ~ 0.061). 
Our present result in the low doping range also is well 
consistent with this recent experimental observation^. 



S = 0.06 (solid line), S = 0.09 (dashed line), S = 0.12 
(dash-dotted line), and S = 0.15 (dotted line) in Fig. [5]in 
comparison with the corresponding experimental result^ 
of YBa2Cu307_j, (inset). Our result indicates that the 
superfluid density p s {T) decreases with increasing tem- 
perature, and vanishes at the SC transition temperature 
T c . Moreover, the most striking feature of the present 
results is the wide range of linear temperature depen- 
dence at low temperature, extending from close to the SC 
transition temperature T c to down to the temperatures 
T k 4K^ 8K for different doping concentrations. How- 
ever, in correspondence with the nonlinear temperature 
dependence of the magnetic field penetration depth at 
the extremely low temperatures shown in Fig. [3J the su- 
perfluid density p s (T) crosses over to a nonlinear temper- 
ature behavior at the extremely low temperatures (below 
T « 4K^ 8K for different doping concentrations). Our 
these results are also well consistent with the correspond- 
ing experimental result^— for cuprate superconductors. 
The good agreement between the present theoretical re- 
sults and experimental data also shows that the d-wave 
contribution to the superfluid density p s (T) in cuprate 
superconductors is predominant. 




FIG. 4: Doping dependence of the zero-temperature su- 
perfluid density with parameters t/J — 2.5, t' /t — 0.3, 
and J = 1000K. Inset: the corresponding experimental re- 
sult for Yo.8Cao.2Ba2(Cui_ z Zn z )307-<5 (open circles) and 
Tli_j / Pb H Sr2Cai_ a; Y;eCu207 (solid triangles) taken from Ref. 
0. 

The doping dependence of the superfluid density shown 
in Fig. [4] also is strongly temperature dependent. In par- 
ticular, when the temperature T = T c , the kernel of the 
response function K^fa — > 0, 0)|t=t c = and the mag- 
netic field penetration depth A(T C ) = oo as mentioned 
in Subsection IIII Al this leads to the superfluid density 
p s {T c ) = 0, which is consistent with the result of the ef- 
fective superfluid density obtained from Eq. (|23|) . To 
show the superfluid density clearly for all the temper- 
ature T < T c , we plot the superfluid density p s (T) as 
a function of temperature for the doping concentration 



25 




T(K) 

FIG. 5: Temperature dependence of the superfluid density 
for the doping concentration S = 0.06 (solid line), S = 0.09 
(dashed line), S = 0.12 (dash-dotted line), and S = 0.15 
(dotted line) with parameters t/J — 2.5, t'/t — 0.3, and 
J = 1000K. Inset: the corresponding experimental result for 
YBa2Cu307_ !/ taken from Ref. fl5l . 

The essential physics of the nonlinearity in the tem- 
perature dependence of the penetration depth (then the 
superfluid density) in cuprate superconductors at the ex- 
tremely low temperatures, as shown in Fig. [3] (then Fig. 
[5]), in the present t-t'-J model is the same as that in the t- 
J— , and can be attributed to the nonlocal effects induced 
by the gap nodes on the Fermi surface in a pure d-wave 
pairing stated—. A weak external magnetic field acts 
on the SC state of cuprate superconductors as a pertur- 
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bation. Within the linear response theory, one can find 
that the nonlocal relation between the supercurrent and 
the vector potential (j4} in the coordinate space holds due 
to the finite size of Cooper pairs. In particular, in the 
present kinetic energy driven d-wave SC mechanism 22 , 
the size of charge carrier pairs in the clean limit is of the 
order of the coherence length £(k) = hvp / V Ah(k) , where 
vp = h~ 1 d£,k/dk\k F is the charge carrier velocity at the 
Fermi surface, which shows that the size of charge carrier 
pairs is momentum dependent. In general, although the 
external magnetic field decays exponentially on the scale 
of the magnetic field penetration length A(T), any non- 
local contributions to measurable quantities are of the 
order of k~ 2 , where the Ginzburg-Landau parameter k 
is the ratio of the magnetic field penetration depth A and 
the coherence length £. However, for cuprate supercon- 
ductors, because the pairing is d-wave, the charge carrier 
gap vanishes on the gap nodes on the Fermi surface, so 
that the quasiparticle excitations are gapless and there- 
fore affect particularly the physical properties at the ex- 
tremely low temperatures. This gapless quasiparticle ex- 
citation leads a divergence of the coherence length £(k) 
around the gap nodes on the Fermi surface, and then 
the behavior of the temperature dependence of the mag- 
netic field penetration depth (then the superfluid den- 
sity) depends sensitively on the quasiparticle scattering. 
At the extremely low temperatures, the quasiparticlcs 
selectively locate around the gap nodal region, and then 
the major contribution to measurable quantities comes 
from these quasiparticles. In this case, the Ginzburg- 
Landau ratio At(k) around the gap nodal region is no 
longer large enough for the system to belong to the class 
of type-II superconductors, and the condition of the lo- 
cal limit is not satisfied^, which leads to the system in 
the extreme nonlocal limit, and therefore the nonlinear 
behavior in the temperature dependence of the magnetic 
field penetration depth (then superfluid density) is ob- 
served experimentally^^. On the other hand, with in- 
creasing temperature, the quasiparticles around the gap 
nodal region become excited out of the condensate, and 
then the nonlocal effect fades away, where the momen- 
tum dependent coherence length C(k) can be replaced 
approximately with the isotropic one £o = hv-p/nA^. In 
this case, the calculated Ginzburg-Landau parameters 
are k ~ A(0)/Co ~ 166.29, k « 175.55, and k w 156.14 
for the doping concentrations 6 = 0.14, 6 = 0.15, and 
S = 0.18, respectively, and therefore the condition for 
the local limit is satisfied. In particular, these theoreti- 
cal values of the Ginzburg-Landau parameter in differ- 
ent doping concentrations are very close to the range 
Ko w 150 ~ 400 estimated experimentally for different 
families of cuprate superconductors in different doping 
concentrations^^— iSl. As a consequence, our present 
study shows that cuprate superconductors at moderately 
low temperatures turn out to be type-II superconductors, 
where nonlocal effects can be neglected, then the electro- 
dynamics is purely local and the magnetic field decays 
exponentially over a length of the order of a few hun- 



dreds nm. 



IV. CONCLUSIONS 

Within the t-t'-J model, we have discussed the dop- 
ing and temperature dependence of the Meissner effect 
in cuprate superconductors based on the kinetic energy 
driven SC mechanism. Our results show that in the lin- 
ear response approach, the electromagnetic response con- 
sists of two parts as in the conventional superconductors, 
the diamagnetic current, which is the acceleration in the 
magnetic field, and the paramagnetic current, which is 
a perturbation response of the excited quasiparticle and 
exactly cancels the diamagnetic term in the normal state, 
then the Meissner effect is obtained for all the tempera- 
ture T <T C throughout the SC dome. Within this frame- 
work, we have reproduced well all the main features of the 
doping dependence of the local magnetic field profile, the 
magnetic field penetration depth, and the superfluid den- 
sity in terms of the specular reflection model. The local 
magnetic field profile follows an exponential law, while 
the magnetic field penetration depth shows a crossover 
from the linear temperature dependence at low temper- 
atures to a nonlinear one at the extremely low tempera- 
tures. Moreover, in analogy to the domelike shape of the 
doping dependent SC transition temperature, the super- 
fluid density increases with increasing doping in the lower 
doped regime, and reaches a maximum around the criti- 
cal doping S w 0.195, then decreases in the higher doped 
regime. The good agreement between the present theo- 
retical results in the clean limit and experimental data 
for different families of cuprate superconductors also pro- 
vides an important confirmation of the nature of the SC 
phase of cuprate superconductors as a d-wave BCS-like 
SC state within the kinetic energy driven SC mechanism. 

Finally, it should be emphasized that in the present 
study, the only coupling of the electron charge to the 
weak external magnetic field is considered in terms of 
the vector potential A, while the coupling of the elec- 
tron magnetic momentum with the weak external mag- 
netic field in terms of the Zeeman mechanism has been 
dropped. In this case, the above obtained results are only 
valid in the weak external magnetic field limit. How- 
ever, the depairing due to the Pauli spin polarization is 
very important in the presence of a moderate or strong 
external magnetic field, since cuprate superconductors 
are doped Mott insulators with the strong short-range 
antiferromagnetic correlation dominating the entire SC 
phase 4 . In particular, in the the kinetic energy driven 
SC mechanism 2 ^, where the charge carrier-spin interac- 
tion from the kinetic energy term induces a d-wave pair- 
ing state by exchanging spin excitations. Therefore under 
the kinetic energy driven SC mechanism, a moderate or 
strong external magnetic field aligns the spins of the un- 
paired electrons, then the d-wave electron Cooper pairs 
in cuprate superconductors can not take advantage of 
the lower energy offered by a spin-polarized stated. In 
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this case, we 4 ^ have studied the magnetic field depen- 
dence of the superfiuid density and doping dependence 
of the upper critical magnetic field in cuprate supercon- 
ductors for all the temperature T < T c throughout the 
SC dome by considering both couplings of the electron 
charge and electron magnetic momentum with a mod- 
erate and a strong external magnetic field, respectively, 
and the results show that the external magnetic field in- 
ducing an reduction of the low-temperature superfiuid 
density, while the maximal upper critical magnetic field 
occurs around the optimal doping, and then decreases 
in both underdoped and overdoped regimes, in quali- 
tative agreement with the corresponding experimental 
dat a 41 ' 42 . These and the related results will be presented 



elsewhere. 
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Appendix A: Paramagnetic response kernel at zero-temperature and superconducting transition temperature 



In this Appendix, we discuss the paramagnetic part of the response kernel in Eq. (|19p in the long wavelength limit 
at the zero-temperature (T = 0) and SC transition temperature (T = T c ). Firstly, we discuss the case at T = 0. In 
the long wavelength limit, i.e., |q| —> 0, the paramagnetic part of the response kernel (|T9"1) can be evaluated as, 



^ 3 (q^0,0)| T ^ 



-2Z, 



hF' 



4e^ 



N 



sin 2 k y [x\t — 2x2*' cos k x 



j3e 



I) 2 



(Al) 



T^0 



where N — N x N y , N is the number of sites on a square lattice, while N x and N y are corresponding numbers of sites 
in the x and y directions, respectively. In the d-wave pairing state, the characteristic feature is that the energy gap 
vanishes A^k)]^^^ i = Ab(cosk x — cosk y ) /2\ij i ^—i] (: — along the diagonal directions in the Brillouin zone. In 
this case, the paramagnetic part of the response kernel in Eq. (jAlj) can be rewritten as, 



A^)(q^0,0)| T ^o = -2^h F fr 



2 4 e 2 1 



— ^ sm2 k v [Xit ~ 2x2*' cos k x 
k(|fc*l#|fcy|) 

tt sin 2 ky[xit - 2x 2 t' cos k. 



fie 



( e /^hk + 1)2 



{e Kk v + l) 2 



(A2) 



T->0 



where = Z\ l -p(ZtxiCosk y — Zt'xi cos 2 k y — The first term of the right-hand side is equal to zero since the 
existence of the gap, while the second term of the right-hand side can be evaluated explicitly as, 



2Z 2 ^± 
2ZhF H 2 N x 



2Z hF 4e 2 1 
Z ft 2 ~Wx 

2Z hF 4e 2 1 
Z H 2 ~N~ T 



n dk 
.„ 2tt 



y sin 2 k y [x\t — 2x2^' cos k y ] 2 



dk 1 
— ^[xitcosk y - 2x 2 i' cos{2k y )]— ? 

-7T 2-K e Kky _|_ 1 

dk 

-^-[Xit cos k y - 2x 2 i' cos(2k y )}0(£ ky 



J T->0 



T^0 



(A3) 



which is equal to zero in the thermodynamic limit N — > oo (then N x — > oo and N y — > oo), where the step function 
6{x) = 1 for x < and 9(x) — for x > 0. 

Now we turn to discuss the paramagnetic part of the response kernel (|19p in the long wavelength limit at T = T c 
(p c — T^ 1 ). In this case, the energy gap Ahz(k)|T=T c = 0, and the paramagnetic part of the response kernel (fH?|) can 
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be evaluated explicitly as, 

tf#(q->0,0) = _2Z h 2 F ^l^sin 2 fe y [ Xl i-2 X2 t'cosM 2 ' " 



h 2 N ^ 1 ( e P°& + l) 

k 



2 



a 4e 2 r dfc, dk v . 2 , a /3 c e^ 

~ 2Z hF^T J -Tfr J -r^-sm k y [xit-2 X 2t cosk x \ 



.2 l-TT 



„ 4e^ r dk x r dk Vr , „ , , , . 1 

= hF 7^ J_ v ^-[Xitco8k y -2 X2 t cosfc, cos fc y ] e/3e& + 1 

4e 2 1 ^-^ - 4e 2 1 

= Z hF — — ^Jxiicos/c^ _ 2x 2 i'cosfc x cosfc y ]n F (^k) = -^-[Xi^i* ~ 2 X 2 <M'] = -^5", (A4) 



which exactly cancels the diamagnetic part of the response kernel in Eq. (|9]), then the Meissner effect is obtained for 
all T <T C throughout the SC dome. | 
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